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Abstract 

We report a preliminary measurement of the D^-D " mixing parameter ycp and the CP-violating 
parameter using the decay D*~^ D^tt^ followed by K~tt~^ and K~^K~ . The 

results are obtained from a 158 fb""*^ data sample collected near the T(45') resonance with the Belle 
detector at the KEKB asymmetric energy e^e~ collider. 
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INTRODUCTION 



The rate of D^-D^ mixing, which has not yet been observed, would provide an impor- 
tant window on new physics. The contribution of box diagrams is very small due to GIM 
suppression, and the mixing rate is believed to be dominated by long-distance processes. 
These processes are themselves suppressed by S'?7(3)-flavor symmetry. The mixing rate is 
usually expressed in terms of the parameters x = Am/T and y = Ar/(2r), where Am and 
Ar are the mass- and width-differences of thephysical states, and F is the average decay 
width. Numerous predictions for x and y exist a recent calculation gives 10~^ j2[. 

Non-Standard-Model (SM) processes are expected to enhance x relative to y, and measuring 
X ^ y would be a strong indication of new physics 0, 1^ . 

In this paper we report a measurement of the parameter 

_ T{K~'n+) 

which is equal to y in the limit where CP is conserved. If CP is violated, the widths of 
and decays into the CP eigenstate K^K~ differ. These widths can be expressed as 

V{D^ ^ K+K-) = T [1 + Rm{y cos (j)-x sin (j))] 
t(D^ ^ K+K-) = t[i + R-\y cos (p + x sin 4 

where and are CP- violating parameters defined in j^]. CP violation in mixing results in 
Rm 7^ 1, and CP violation due to interference between mixed and unmixed decay amplitudes 
results in 7^ 0. We investigate whether CP violation is present by measuring the quantity 

The SM predicts A-p ^ 10~^ [5]; a large value would indicate physics beyond the SM. 

Our measurement of ycp is obtained by measuring the lifetime difference between — *■ 
K~TT'^ and — > K^K~ decays. (Throughout this paper, charge-conjugate modes are 
included unless noted otherwise.) Our measurement of is obtained by measuring the 
lifetime difference between _D° K^K" and — > K^K~ decays. The analysis is based 
on 158 fb~^ of data, and all results are preliminary. To identify the flavor of the and also 
reduce background, we require that all D°'s originate from D*^ — > D^n^ decays. 



EVENT SELECTION 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three- 
layer silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), an array of 
aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-fiight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals 
(ECL) located inside a super- conducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to detect mesons and to 
identify muons (KLM). The detector is described in detail elsewhere 

We select well-measured charged tracks that have impact parameters less than 1.0 cm in 
the radial direction and less than 2.0 cm along the beam direction [z) with respect to the 
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interaction point (IP). The daughter tracks are required to have at least two hits in each 
of the r-(j) and z-measuring SVD layers. This requirement is not applied to the slow pion 
from D*~^ — >■ D^7r~^ decays. 

For charged particle identification, information from the CDC, TOF, and ACC subsys- 
tems is combined to form an overall hadron likelihood: 

C{h) = C'-'^'^ih) X /:CDC(;^) ^ ^TOF^^) ^ 

where h stands for vr or K. Charged particles are identified as pions or kaons using the 
likelihood ratios: 

V(K) = ^ Vin) = ^ = l-V(K) 

Kaon candidates are required to have V{K) > 0.6. This requirement is 88% efficient and 
has a pion misidentification probability of 8%. Pion candidates are required to satisfy 
V{K) < 0.9. 

Both continuum and T(4S') data are used in this analysis. Each event is required to 
contain at least one D*~^ meson decaying to D^n^. To remove D*~^^s resulting from B meson 
decays, the D*^ is required to have a momentum in the e^e~ center-of-mass frame (CMS) 
greater than 2.5 GeV/c. meson candidates are reconstructed from K~tt^ and 

DO ^ K+K- decays. All D° candidates must satisfy \M{D^){it-M{K-7i+ /K+K-)\ < 2a&t, 
where M(D°)fit and (Tfit are the mean and standard deviation of a Gaussian fit to the mass 
peak in data. These parameters are calculated independently for each decay mode in 
order to help keep the K~7r~^ and K^K~ vertex resolution functions unbiased by the mass 
requirement. From the fits, (Jfit is found to be 5.1 MeV/c^ for K^tt^ and 4.5 MeV/c^ 

for — > K^K~ (see Table |l}. The momentum of the slow from D*^ —>■ D^tc^ decays is 
obtained by refitting using the beam position constraint, i.e., the track must project back to 
the IP region. The invariant mass of the D*'^ candidates after slow 7r+ refitting is required 
to satisfy \M{D\+) - M {R- 11+ / K+ R-) + M^^^ - MlV^\ < 1.0 MeV/c^. 



VERTEX FIT AND DECAY TIME RECONSTRUCTION 



The proper time between the production and decay of the candidates is determined 
by projecting the momentum P and the flight length £ onto the xy plane transverse to the 
beam axis: 

^xy ■ Pxy ( T^no \ 

= ^ ^ (— ) ■ 

The momentum P is the vector sum of the R and vr momenta, and £ is the displacement 
vector between the production and decay vertices. Information from the z projections is not 
included due to the large uncertainty in the z component of the IP. 

To improve the spatial resolution, the vertex fit of — > R~7f^{R^R^) decays is done 
in three steps. First, the decay point is fit using only information from the daughter 
tracks. Then, the production vertex is found by fitting the reconstructed D° momentum 
vector with the IP region, whose profile is obtained from a separate study of the full data 
set. Finally, the decay vertex is corrected using information from the production 
vertex fit. No mass constraint is applied to the D° candidates during the vertex fit to avoid 
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mode-dependent bias that could result from a correlation between the invariant mass and 
a global momentum correction procedure. The error in the position of the D° decay vertex 
obtained from the fit is required to be less than 150 fim in each of the x and y directions. 
Fig. [T] illustrates the high efficiency of this cut. This procedure produces typical D° decay 
vertex resolutions of 45 fim in both x and y. 
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FIG. 1: Distributions of decay vertex errors ax (left plot) and ay (right plot) returned by the 
vertex fit after the IP position correction. The horizontal scale is in cm. Vertical lines indicate the 
position of the fit quality cut: ax^y < 150 fim. 
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FIG. 2: Invariant mass spectra for — > K~-k~^ candidates (left plot) and K^K~ candidates 

(right plot). The horizontal scale is in GeV/c^. The bin size is 0.5 MeV/c^ and the solid curves 
show the fit results. 

The invariant mass spectra of the selected K^t^^ and K^K~ candidates 

are shown in Fig. |21 The distributions are fit with a sum of two Gaussians representing the 
signal and a first-order polynomial representing the background. The results of the fit are 
summarized in Table B 



RESOLUTION FUNCTION STUDY 



The proper time distribution is represented by an exponential function with time constant 
Tj^ convolved with the detector vertex resolution. The properties of the detector resolution 
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TABLE I: Invariant mass spectra fit results. 



Mode 


# events 


main Gaussian a 


main Gaussian fraction 


Signal purity 






(MeV/c2) 


(%) 


(%) 




448000 


5.1 


78.4 


99.1 


K+K- 


36480 


4.5 


74.0 


97.6 
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FIG. 3: Proper time resolution functions for —>■ K~'ir~^ (left plot) and — > K~^K~ (right plot) 
obtained from MC simulation. The horizontal scale is in fs. The bin size of the data points is 40 fs 
and the solid curves show the results of the fit. (The fit is performed using a bin size of 10 fs.) 

function are studied using Monte Carlo (MC) simulation of the signal decay modes. Ap- 
proximately 7 ■ 10^ — > K'TT^ decays and 1 ■ 10^ — >■ K^K~ decays were simulated, 
corresponding to two times the data sample and three times the data sample, respectively. 
Fig. 121 shows the difference between the generated and reconstructed decay times of the 
mesons. The distributions are fit with a sum of five Gaussians constrained to a common 
central value. The results of the fit are summarized in Table In] the confidence level (CL) 
is 22.2%. Due to similar kinematics, the K^K~ decay mode is fit with the same resolution 
function as that used for K^-k^ . However, to account for small differences between K^K~ 
and i^~7r"'", the width of each Gaussian component in the K^K~ fit is multiplied by a 
common scale factor a. Fitting first the K~tt^ distribution and then the K^K~ distribution, 
we obtain a = 1.043±0.004. The CL for the K^K~ fit is 71.3%. The common central values 
of the Gaussians are found to be Xq = —1.51 ± 0.22 fs for K~n~^ and Xq = —1.57 ± 0.63 fs 
for K+K-. 

In the lifetime difference fits for both MC and data, the relative fractions of the Gaussian 
components of the resolution function are fixed to the values determined in the above study. 
The central value of the K^K~ resolution function is fixed to that of the K'lr^ resolution 
function, which is consistent with the results above. The possible bias due to this assumption 
is included in the systematic error. 

MC STUDY OF LIFETIME DIFFERENCE FIT 

The lifetime difference is obtained from a simultaneous, binned, maximum likelihood 
fit to the measured proper lifetime distributions of — * K^K'^ and K^ti^ decays. 
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TABLE II: D K vr^ proper time resolution function fit results. 



Fit parameter 


Fraction (%) 


Value (fs) 


0-1 


26.1 


95.1 ± 1.3 


0-2 


50.4 


177.0 ±2.2 


0-3 


19.8 


328.7 ±7.4 


£74 


3.1 


675.7 ±24.9 


0-5 


0.6 


2199 ± 95 



The consistency and quality of the fitting procedure is tested with MC events. The lifetime 
distributions of — *■ Kt: /KK decays were fit with an exponential function convolved with 
a resolution function of the form described above. The CL of the fit is 61%. Fig. E] and 
Table UTTl show the results of the fit and a comparison to the parameters obtained from the 
resolution function study. The parameters are seen to agree to within the quoted errors. The 
resolution function scale factor a determined from the fit is 1.042 ±0.007, also in agreement 
with the result from the resolution function study. The lifetime obtained is 411.9 ±0.8 fs, 
which is consistent with the generated value of 412.6 fs. The lifetime difference obtained is 
0.0 ± 1.7 fs, which is consistent with the generated value of exactly fs. 




-4000 -2000 2000 4000 -4000 -2000 2000 4000 

FIG. 4: Simultaneous lifetime difference fit for K~-k~^ (left plot) and — > K~^K~ (right 

plot) performed on MC-generated events. The horizontal scale is in fs. The bin size of the data 
points is 40 fs and the solid curves show the results of the fit. (The fit was done using a bin size 
of 10 fs.) 



LIFETIME DIFFERENCE FIT IN DATA 

To fit the lifetime difference in data, one must include small contributions from back- 
ground processes. To include these, the proper time distribution of background is obtained 
from events in the mass sidebands. The sidebands chosen begin ±25 MeV/c^ from m^o 
and cover the same total range as the signal region (4crfit)- 

The background proper time distributions are parameterized as the sum of an exponential 
convolved with a Gaussian and a Breit-Wigner function. The resulting spectra are shown in 
Fig. The systematic uncertainties due to the background fitting procedure are evaluated 
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TABLE III: Results of the simultaneous lifetime fit (MC) vs. standalone resolution function fit. 



Fit parameter 


Resolution Function (fs) 


Lifetime Fit (fs) 


0-1 


95.1 ± 1.3 


94.4 ± 1.7 




1 < ( .U it Z.Z 


1 70 n 1 1 
i < y.u ± i.z 




328.7 ± 7.4 


328.2 ±2.2 


f^i 


675.7 ±24.9 


(364.4 ± 8.5 


0-5 


2199 ± 95 


2225 ± 70 


^0 


-1.51 ±0.22 


-0.95 ± 0.54 


o 


1.013 ±0.001 


1.0-12 ±0.007 
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FIG. 5: Background proper time distributions for K tt'^ (left plot) and K^K (right 

plot). The horizontal scale is in fs. The bin size is 40 fs and the solid curves show the results of 
the fit. 

by varying the fitting parameters. Since the background contribution is small compared to 
the signal, the effect of the background magnitude and shape uncertainty is small (approx- 
imately 0.1% in Ucp)- It was checked with MC simulation that the small fraction of signal 
events remaining in the mass sidebands does not affect the ycp measurement above the 
quoted systematic error on background. The decay time function obtained for background 
is included (with normalization fixed) in the lifetime fit to the signal region. 



TABLE IV: Results of the simultaneous lifetime fit (data) vs. standalone resolution function fit. 



Fit parameter 


Resolution Function (fs) 


Lifetime Fit (fs) 


0-1 


95.1 ±1.3 


112.3 ±2.6 


0-2 


177.0 ± 2.2 


198.2 ± 1.8 


0-3 


328.7 ±7.4 


378.0 ± 3.2 


(74 


675.7 ±24.9 


864.2 ± 16.4 




2199 ± 95 


3197 ± 309 


Xo 


-1.51 ±0.22 


8.6 ±0.8 


a 


1.043 ± 0.004 


1.056 ±0.012 
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FIG. 6: Simultaneous lifetime difference fit to the data for K'tt+ (left plot) and 

K^K~ (right plot). The horizontal scale is in fs. The bin size is 40 fs and the solid curves show 
the results of the fit. The dashed curves show the background component included in the fit. 

Fig. ini and Table IIVI show the results of the lifetime difference fit to the data. The CL 
of the fit is 94%. The fitted K-7i+ lifetime is 412.6 ± 1.1 fs, and the fitted lifetime 

difference between K^-k^ and K^K~ is 4.76 ± 2.85 fs. Together these values 

give ycp = 1.15 ±0.69%. 

The parameters obtained by the fit for the resolution function (see Table IIV|1 differ from 
those obtained from the MC study by more than the quoted errors. This is nominally due 
to vertex detector misalignment, IP mismeasurement, material unaccounted for in the MC 
simulation, and momentum dependence. These effects are included when evaluating the 
overall systematic error on y^p (Table IVj) . 



TABLE V: Systematic errors of the ycp measurement. 



Source 


Value (%) 


Xo(-fC7r) = Xq{KK) assumption 


0.15 


a difference on data and MC 


0.20 


Bin size 


0.10 


mass window 


0.20 


Background 


0.10 


Selection requirements 


0.15 


Total 


0.38 



CPV ASYMMETRY FIT 

The CP-violating parameter in D''' K^K~ decays is obtained from fits to the data 
using the charge of the slow pion from D*^ — >■ D^tt^ or D^7t~ to distinguish decays 
from D ° decays. The and D ^ subsamples are fit separately using the resolution function 
parameters obtained from the full statistics fit. Background is treated as described in the 
previous section. The same procedure is performed for K~tt^ decays and also for 
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MC events in order to estimate the systematic uncertainty of the fit. All results are listed 
in Table IVll The systematic uncertainty of the measurement is taken to be ±0.30%. 



TABLE VI: Results of on data and MC. 



decay mode 


data/MC 


^(,=+1) - ^(,=-1) (fe) 


(%) 


Ktt 


MC 


1.2 ± 1.1 


0.15 ±0.13 


KK 


MC 


2.6 ±3.1 


0.31 ±0.38 


Ktt 


data 


-0.76 ± 1.55 


-0.09 ±0.19 


KK 


data 


-1.67 ± 5.18 


-0.20 ±0.63 



SUMMARY 

In conclusion, we report a preliminary measurement of the D^-D^ mixing parameter ycp 
and the CP-violating parameter obtained by measuring lifetime differences in 

and K+K' decays. The flavor of the D° or D° is identified via D*+ D°7i+ 

decays. The results obtained are ycp = (1.15 ± 0.69 ± 0.38)% and Aj. = (-0.20 ± 0.63 ± 
0.30)%. 
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